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Abstract 


GRB 200612A could be classified as an ultralong gamma-ray burst due to its prompt emission lasting up to 
~1020s and the true timescale of the central engine activity tours, > 4 X 10*s. The late X-ray light curve with a 
decay index of a = 7.53 is steeper than the steepest possible decay from an external shock model. We propose that 
this X-ray afterglow can be driven by dipolar radiation from the magnetar spindown during its early stage, while 
the magnetar collapsed into the black hole before its spindown, resulting in a very steep decay of the late X-ray 
light curve. The optical data show that the light curve is still rising after 1.1 ks, suggesting a late onset. We show 
that GRB 200612A’s optical afterglow light curve is fitted with the forward shock model by Gaussian structured 
off-axis jet. This is a special case among GRBs, as it may be an ultralong gamma-ray burst powered by a magnetar 


in an off-axis observation scenario. 


Key words: (stars:) gamma-ray burst: individual (GRB 200612A) — (stars:) gamma-ray burst: general — 
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1. Introduction 


As the most energetic events in astronomy, gamma-ray 
bursts (GRBs) have been extensively studied. The GRBs 
whose prompt emission duration is less than 2 s are classified as 
short-GRBs, while those lasting for more than 2s are called 
long-GRBs. In the histogram of duration, short GRBs peak at 
0.2-0.3 s, while long GRBs peak at 20-30 s (Kouveliotou et al. 
1993) and ranging from ~2 s up to several hundred seconds in 
the observer frame. 

Some GRBs have been reported with an extraordinary 
duration of 10* s, known as ultralong GRBs. It has been 
debated whether these belong to a new distinct class of events 
or if they are similar to long GRBs. The former hypothesis 
suggests that the progenitors of these ultralong GRBs may be 
blue supergiants (Gendre et al. 2013; Nakauchi et al. 2013; 
Levan et al. 2014; Huang et al. 2020), while the latter 
hypothesis proposes that ultralong GRBs share the same 
progenitors with long gamma-ray bursts, representing the tail of 
the distribution of long GRBs (Virgili et al. 2013; Evans et al. 
2014). Zhang et al. (2014) proposed defining ultralong GRBs 
based on the true timescale of the central engine activity (fpurst)- 
The fpurst distribution of their “good” sample shows an apparent 
bimodal distribution with a narrow, significant peak at ~355 s, 
and a wider, less significant peak at ~2.8 x 104s, respectively. 
Gao et al. (2017) obtained similar results as a narrow, 
significant peak at 173.8 s and a wider, less significant peak at 
~2.75 x 10*s , respectively. As an example of an ultralong 
case, GRB 101225 shows that plateau and extremely sharp flux 


drop could present in X-ray afterglow, it may be driven by a 
newborn magnetar (Lii et al. 2018; Zou et al. 2019). Both the 
normal long GRB model and the blue supergiant model will 
launch an ultra-relativistic jet. Usually, we assume that 
observations are occurring on-axis, but some off-axis GRBs 
have been observed, such as GRB 170817A (Abbott et al. 
2017). An observer at 0, < 6;, where 0, is the viewing angle, 0; 
is the jet half-opening angle, should see a light curve very 
similar to that for an on-axis observer, an observer at 0, > 0; 
will see a rising light curve at early times, peaking when the jet 
Lorentz factor is ~1/6,, and approaching that is seen by an on- 
axis observer at later times (Woods & Loeb 1999; Granot et al. 
2002). 

GRB 200612A is a special case that may be an ultralong 
gamma-ray burst driven by amagnetar spinning down in an 
off-axis scenario. The structure of the paper is as follows. The 
observations are presented in Section 2. In Section 3, we 
perform the analyses of afterglow. Then, in Section 4, we 
model the light curves. Our conclusions and discussion are 
given in Section 5. 


2. Observations 


GRB 200612A triggered the Swift Burst Alert Telescope 
(BAT) on 2020 June 12 at 04:39:37 UT (denoted as Tọ in this 
paper). The mask-weighted light curve shows a complex 
structure that starts at ~To — 20s and lasts until ~7) + 200s. 
Too (15-350 keV) estimated by BAT is 
Too, BAT = 171.48 + 19.74s (Markwardt et al. 2020; Sonbas 
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Figure 1. The prompt emission light curves of GRB 200612A with 4.096 s 
binning. Observations from the Fermi/GBM detectors and Swift/BAT 
detectors are represented by red solid and blue solid lines, respectively. 


et al. 2020). At 04:49:38.68 UT (To + 10 minutes), the Fermi 
Gamma-Ray Burst Monitor (GBM) triggered and located 
GRB 200612A. Fermi GBM Team report they detected this 
burst and provide To) Gam ~ 88s (50-300 keV) at GCN 
circular. In Fermi GBM Burst Catalog,' they provide 
Too,cBM = 984+ 65s. The light curve of Fermi/GBM and 
Swift/BAT is presented in Figure 1. GRB 200612A was also 
detected by Konus—Wind (Svinkin et al. 2020). They observed 
that the burst light curve shows four emission episodes in the 
interval from ~Tọ to ~Tọo + 1020 s. They classified this burst as 
an ultralong GRB. As observed by Konus—Wind, the burst had 
a fluence of 4.081};94 x 1075 erg cm~?, 

The Swift X-Ray Telescope (XRT) began observing the field 
at 04:41:04.7 UT (To + 87.7 s). The enhanced Swift/XRT 
position is R.A., decl. (J2000) = 201048249, —45°20/07! 8 
with an uncertainty of 1”5. The Swift Ultraviolet /Optical 
Telescope (UVOT) took a finding chart exposure of 150s with 
the white filter starting at Tọ +95 s. However, no credible 
afterglow candidate has been found in their initial data 
products. 

MASTER-OAFA robotic telescope located in Argentina 
(OAFA observatory of San Juan National University) was 
pointed to GRB200612A at 2020 June 12 04:47:52 UT 
(To +495 s). They found an optical transient within Swift 
error-box (R.A. = 302.679, decl. = —45.3569, r=0.05) with 


1 https: //heasarc.gsfc.nasa.gov /cgi-bin/W3Browse/w3hdprods.pl? 
files=P&Target=heasarc%Ffermigbrst%7C%7C%7C %5F%5Frow=586% 
7C%7C&Coordinates=Equatorial& Equinox=2000 
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Figure 2. The R-band image of GRB 200612A taken on 2020 June 14, with 
four frames combined. Echo frame exposure 300 s was taken from 00:45:34 to 
01:01:58 UTC, i.e., 1.84 days after the trigger of the GRB. 


17.9 Mag at To + 1050s (Lipunov et al. 2020) and 17.5 Mag at 
To+ 1172s (Tiurina 2020). The 60cm BOOTES-3/YA 
robotic telescope at NIWA Lauder in Otago (New Zealand) 
began observing the position reported by Swift/XRT at 
09:51:32 UT(~5.2 hr after trigger). In the co-added image 
(20 exposures of 60 s each in the clear filter), no source is found 
at the position, and the upper limit is 20.3 mag (Hu et al. 2020). 
We observed the field of GRB 200612A with the LCO 1m 
Sinistro instrument starting at 2020 June 14 00:45:34 UT 
(To + 1.84 days), four images in the R-band were obtained 
(exposures 300s each). A weak source was found within the 
Swift error box. Figure 2 shows the R—-band image of the field 
of GRB 200612A taken at Ty + 1.84 days. The data reduction 
was carried out following the standard routine in the IRAF 
package (Tody 1986). The photometric result is 
R= 22.90 + 0.35, not corrected for the Galactic extinction. 


3. Data Analysis 
We fit the X-ray light curve with a single power-law function 
F = Fot~°, (1) 


and a broken power-law function 


t way i war ]-1/w 
F= F P F P3 A (2) 
b b 


where F is the flux at the break time (tp), a, and a» are the 
afterglow flux decay indices before and after tp, respectively, 
and w is a smoothness parameter which represents the 
sharpness of break. Five phases comprise this X-ray light curve 


(Figure 3). Phases I and II form a flare with a; = =1,56"3 1, 
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Figure 3. Broadband light curve of GRB 200612A, including y-ray and X-ray 
flux and optical flux densities. The red solid line represents the best-fitting 


result for the X-ray afterglow using a power-law or broken power-law model. 
The blue point is the X-ray afterglow of GRB 101225A. 


az = 3.061%-19 and peak at t, = 14513 s, Phase III rises from 
~300 s, then Swift satellite suffers from Earth occultation, so 
there are no data between 380 and 3700 s. Phases IV and V are 
fit to broken power-law function with ag = 1.1807. 
as = 7.53*°") and a break time at tp = 3.09977 x 104s. 
The steep decay in Phase V is inconsistent with an external 
shock model. The steepest possible decay from an external 
shock occurs after a jet break. The post-jet-break decay slope 
can be as steep as p, where p is the electron spectral index 
(Rhoads 1999; Sari et al. 1999). Afterglow modeling suggests 
that p is typically between 1.5 and 3; however, as = ie Ut 
obviously deviates from this range. According to the definition 
Of fburst proposed by Zhang et al. (2014) and Gao et al. (2017), 
we obtained fpurst > 4.01 x 10* s of this GRB, aligning with the 
definition of ultralong GRBs. 


4. Model and Result 


The X-ray light curve of GRB 200612A breaks to a steep 
decay at t,, which is far too steep for a jet break. A rapid drop 
of the X-ray flux with a decay index steeper than —3 followed 
by the plateau in some cases which is considered as evidence of 
a magnetar collapsing into a black hole (Troja et al. 2007). We 
propose that a newborn magnetar powered this X-ray emission, 
and it collapsed to form a black hole at t,, resulting in the decay 
of Phase V is extremely steep. The newborn magnetar 
spindown by losing its rotational energy through two channels 
(electromagnetic dipolar radiation and gravitational wave 
radiation Shapiro & Teukolsky 1983), the electromagnetic 
emission luminosity would evolve with time as 
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L œx (1+ t/t)“, where a= 1 when gravitational-wave emis- 
sion dominates the spindown losses, and a@=2 when 
electromagnetic losses dominant instead (Dai & Lu 1998; 
Zhang & Mészáros 2001). Usually, the stage where gravita- 
tional-wave emission dominates the spindown energy lost 
should end in early time. There is more likely electromagnetic 
dipolar radiation dominating before the light curve breaks at tp. 
Therefore we describe this X-ray light curve as follows: 


2 

t 

= t< tp 
ri + 5] 


—2 —a 
tp t 
L 1 + — — , t> tb 
wo 7 Cbs 5) e 


where Lyp is the characteristic spindown luminosity, and T is 
the characteristic spindown timescale. Since this burst lacks a 
measured redshift, we assumed its redshift to be z = 2.27, based 
on the average redshift obtained from the long GRB sample 
(Campisi et al. 2010). 

There is sparse data in the optical band; however, we can still 
infer that the onset is later than 1.1 ks. The onset time 
distribution ranges in the GRBs population from 30 to 3000s 
(Liang et al. 2013). The off-axis jet model may explain the 
observed later onset afterglow (Granot et al. 2002; Sato et al. 
2021). The jet with a large Lorentz factor has a strong 
relativistic beaming effect, making the afterglow emission very 
dim for the off-axis observer. As the jet decelerates, the 
beaming effect weakens, leading to a rising in afterglow light 
curves. We attribute GRB 200612A optical emissions to be 
produced by external forward shock in an axisymmetric 
Gaussian structured off-axis jet case. The external forward 
shock model (Sari et al. 1998; Huang et al. 1999) requires a few 
parameters to calculate the afterglow in multiband. These 
parameters include the isotropic-equivalent kinetic energy 
Ex iso, the initial bulk Lorentz factor I'o, the ambient interstellar 
matter (ISM) which is assumed to be uniform with the number 
density mo, the fractions of shock energy that go to electrons 
and magnetic fields (e. and €g), the electron spectral index p, 
the jet opening half-angle 0; and viewing angles 6,. We 
consider the Gaussian structured jet as: 


Loo 7 


L(t) = > 6) 


02 
EXK.iso(0) = Ek iso,0 a| -2) (4) 
02 
Miso 0) = Miso exp] ——> E 5 
(0) „0 exp 7E (5) 


where Ex iso(9) and Miso(0) represent the isotropic-equivalent 
kinetic energy and isotropic-equivalent masses of the jet at 0, 
respectively. The initial Lorentz factor at 0 is set as 
To(O) = Ex.iso(@)/Miso(@)c? + 1. 8 is the polar angle with 
respect to the jet axis, 0e is characteristic angle for the jet. 
Ex iso(9) and Miso(0) are zero when 0 is greater than the 
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Figure 4. The posterior distributions of parameters for the best-fit light curves of GRB 200612A. 


truncation angle 6,,. We performed a Markov Chain Monte 
Carlo exploration of the multidimensional parameter space 
using emcee (Foreman-Mackey et al. 2013). The one- and two- 
dimensional marginalized views of the posterior distribution for 
the Gaussian jet model fits show in Figure 4. The parameter 
estimation priors and marginalized posteriors are reported in 
Table 1, the optimal fitting results are shown in Figure 5. Both 
X-ray and optical data comply with the model lines. 


5. Conclusions and Discussion 


GRB 200612A prompt emission lasting ~1020 s and the true 
timescale of the central engine activity fpurst > 4 x 104 s, could 


be classified into an ultralong gamma-ray burst. The X-ray light 
curve drops rapidly with a decay index of ~ —7 post break time 
tp ~ 3 x 10*s. This is the signal of a center engine turns off. 
The newly born magnetar dipolar radiation could power this 
X-ray before ¢,, then collapsing into a black hole, making the 
light curve decay sharply. 

This burst’s redshift has not been directly measured. We 
assume its redshift to be z=2.27. We estimated the 
characteristic spindown luminosity to be Lyp ~ 1.7 x 1048 
erg st. while the characteristic spindown timescale is 
estimated to be 7 ~ 6.4 x 10° s. The collapse into a black hole 
is projected to occur at t, ~ 2.8 x 10*s. With the fluence of 
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Table 1 
Parameter Estimation Priors and Marginalized Posteriors for the Afterglow of 
GRB 200612A 


Parameter Unit Prior Form Bounds Posterior 
logiy Lup erg s | uniform [30, 60] 48.22030 
logio T s uniform [3, 5] 3.81181 
logo tp s uniform [3, 5] 4.457033 
m = uniform [1, 10] 6.614054 
logo ĒK,iso,0 erg uniform [47, 59] 54.114939 
1081o Miso,0 g uniform [25, 35] 31.29+0-15 
logyo €e = uniform [—8, 0] — 1.725018 
logig €B X uniform [-8, 0] —1.40104] 
0. degree uniform [0.1, 30] 0.791639 
Oy degree uniform [0.1, 60] 5.73413 
0, degree uniform [0, 30] 1.11+8:80 
logion cm? uniform [-5, 5] 1.11787 
p A uniform [2, 3] 2.351005 
—- XRT spindown 
10° F } =-= XRTFS 10° 
— XRT total 
-- RFS > 
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Figure 5. Modeling the optical and X-ray light curves with magnetar spindown 
adds external forward shock model in off-axis case. The dashed—dotted and 
dashed lines in blue are dipole radiation of magnetar spindown and off-axis 
external forward shock model at X-ray (0.3-10 keV), respectively, the blue 
solid line is the total flux calculated from our model. The red solid line is flux 
density calculated from the off-axis external forward shock model. 


4.081494 x 1075 erg cm~? (Svinkin et al. 2020), the isotropic- 
equivalent luminosity will be L, iso = 4.8°97 x 10° erg s~!. 
The isotropic X-ray energy released during the magnetic dipole 
radiation epoch is Ex iso = Lupt/(1 + z) = 1.4483 x 105 
erg. One can define the radiative efficiency for a magnetar to 
convert its spindown energy to radiation as 
x = Ex iso/(Ex,iso + Ex.iso), then obtain hk = 0.01 10.007" 
Furthermore, we estimate the strength of the dipolar magnetic 
field as By 15 = 2.05(l4sR¢ *Lyp/joT3 ) G and the initial spin 
period as P),-3 = 1.42(jf7Lyyp/foT3 ^) 8, where Q, is Q/10" 
in the cgs units. Taking the moment of inertia 7 = 10” g cm”, 
and R=10°cm, we derive B, = 0.7863] x 10'°G and 
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Figure 6. The isotropic luminosity vs. the duration time of different classes of 
transients. The red star is GRB 200612A with assumed redshift z = 2.27. 
The magenta dots are ultralong GRBs taken from Huang et al. (2020), they 
collect these from (Nicastro et al. 2004; Campana et al. 2006; Starling 
et al. 2011; Zhang et al. 2012; Virgili et al. 2013; Zauderer et al. 2013; Evans 
et al. 2014; Cucchiara et al. 2015; Frederiks et al. 2016; Hou et al. 2018), the 
orange dots are ultralong GRBs taken from Levan et al. (2014), and other 
transient data are taken from Levan et al. (2014). 


P) = 1.371038 ms. These values are in good agreement with 
those reported by Lü & Zhang (2014). 

The optical data show that the light curve still rising post 
~1.1 ks, has a late onset, this is in accord with the character of 
the off-axis observer case. We model the optical afterglow 
using an off-axis Gaussian structured jet model and find that 
the viewing angle is 0, ~ 1°1. The truncation angle 6,, ~ 5°7 is 
great than 0,, allowing the prompt emission to be observed. 

As shown in Figure 6, although the duration To of 
GRB 200612A is much smaller than several famous ultralong 
GRBs, for instance GRB 101225A (To) > 7000s; Levan et al. 
2014), GRB 111209 (To) ~ 10000 s; Golenetskii et al. 2011), 
GRB 121027A (To) ~ 6000s; Levan et al. 2014), but it still 
located in same zone of the ultralong GRBs reported at (Huang 
et al. 2020), e.g., the magenta dots in Figure 6, and nearby the 
ultralong GRB 091024A (To9 = 1300s, Virgili et al. 2013). 
GRB 101225A is a famous ultralong GRB with To9 > 7000s 
and fpurst > 1.07 x 10°s (Campana et al. 2011; Thöne et al. 
2011; Nakauchi et al. 2013; Levan et al. 2014; Zhang et al. 
2014; Gao et al. 2017). Its X-ray exhibits a normal power-law 
decay segment (Fart '*'*°°) followed by an extremely 
sharp flux drop (Far °**°4), The possible origin of 
ultralong GRB 101225A may be consistent with the magnetar 
spin-down emission dominated and collapses to a black hole at 
a later time (Lü et al. 2018; Zou et al. 2019). The decay index is 
consistent with Phase IV (ag= ee Ec i) and V 
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(as = 7.53973) of GRB200612A with the error. The 
similarity of the light curves suggests that GRB 200612A and 
GRB 101225A may have the same origin. 

Classifying GRBs by duration is rough. Virgili et al. (2013) 
take the critical duration as 1000s, while Levan et al. (2014) 
suggest that an ultralong burst high-energy emission should be 
lasting to ~10*s. However, the origin of ultralong bursts 
remains a subject of debate. Tidal disruption events (Shcher- 
bakov et al. 2013; MacLeod et al. 2014) and a GRB with a blue 
supergiant progenitor (Nakauchi et al. 2013) are potential 
explanations for these ultralong bursts. However, Virgili et al. 
(2013) propose that ultralong GRBs may represent the tail of 
the duration distribution of the long GRB population. In this 
study, we do not attempt to ascertain the progenitor of 
GRB 200612A. We utilize dipolar radiation from magnetar and 
external shock radiation from an off-axis jet to model the 
afterglow light curve of this burst. 
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